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Abstract. In this paper, we apply an analytical model [V.V. Kulagin et al., Phys. Plasmas 14,113101 
(2007)] to describe the acceleration of an ultra-thin electron layer by a schematic single-cycle laser pulse 
and compare with one-dimensional particle-in-cell (ID-PIC) simulations. This is in the context of creating 
a relativistic mirror for coherent backscattering and supplements two related papers in this EPJD volume. 
The model is shown to reproduce the ID-PIC results almost quantitatively for the short time of a few laser 
periods sufficient for the backscattering of ultra-short probe pulses. 

PACS. 41.75.Jv Laser-driven acceleration, - 52.59.-f Intense particle beams and radiation sources in 
physics of plasmas, - 29.25.-t Particle sources and targets 



1 Introduction 

Irradiation of ultra-thin solid foils by high-contrast laser 
pulses at relativistic intensities may provide a new way 
to create a relativistic mirror for coherent reflection, gen- 
eration of high harmonics and compression of drive and 
probe hght [112131415] . We consider the regime in which 
foil electrons are completely separated from ions _6 . The 
partial reflectivity of such electron layers has been derived 
and compared with one-dimensional particle-in-cell simu- 
lation in [7\ . Coherent reflection depends on the square of 
the layer density and is therefore sensitive to the density 
profile. The analytical model developed by Kulagin et al. 
[T] allows to describe the layer evolution almost quantita- 
tively, at least for a few laser periods after expulsion from 
the foil. This time interval is sufficient to reflect few-cycle 
probe pulses, to compress them to atto- and zepto-second 
duration, and to shift their spectra to the VUV- and X-ray 
regime. Here we apply the Kulagin model to a reference 
case used in the previous publications in order to develop 
and explore this method. 



2 Analytical Model 

Different from [6], we consider here the electromagnetic 
field E = El -I- Eg and B = Bl -I- Bg including the self- 
fields Eg and Bg of the electron layer in addition to the 
external plane laser field, propagating in x-direction and 
having linear polarization with the components Ej^y and 
Blz- Throughout this paper, we use natural units, i.e. 
time and space coordinates are normalized according to 



t' = uj^t, x' = k^x, where ujl are laser frequency and 
wave number, fields E' = eE / {nicLOL), B' = eB / {mcujL}, 
velocities [3 = f/c, and momenta p' — p/mc, where e is 
charge unit, m electron mass, and c velocity of light. In the 
following, these normalized quantities are used dropping 
the prime. 

The self-fields consist of the longitudinal electrostatic 
field Esx due to charge separation between electrons and 
ions and the transverse electromagnetic fields Esy and Bgz 
due to induced electron currents. The laser pulse initially 
hits a foil which has uniform electron and ion density = 
rii — TiQ and thickness do- While the ions are taken as 
immobile heavy particles, the electrons move in the fields. 
Their initial position xo in the layer (0 < a;o ^ f^o) is taken 
as Lagrangian coordinate, and the goal is to determine the 
trajectories x(t,xo) and y{t,xo). 

The longitudinal field, felt by electrons of initial posi- 
tion xq when at position x, is then given by 



Esx{x,xo) 



J A^(a;- 
\N{do 



- a^o) 



(xo < a; < do) 
{x > do) 



(1) 



where 

N ^ no/ricrit ^ (2) 

is the layer density normalized to the critical density ricriu 
it can be conveniently expressed by the squared plasma 
frequency oj^ = Ane'^no/m divided by Equation ^ 
describes the increase of Egx linear in x for electrons still 
inside the ion volume and the constant E^x after leaving 
the ion layer. It is valid as long as electron sub-layers with 
different xq keep their relative order. 

Next we consider the electromagnetic fields due to the 
currents induced in the electron layer by the driving laser. 
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An electron, that originated from xq and has position 
x(t,xo) and velocities (3x(t,xo), (3y{t,XQ) at time t, ex- 
periences approximately 



Esy (t, Xo) 



Bsy{t,Xo) 



Nf3y 


Xo do - Xo 


2 

NPy 


1-/3, 1 + 

Xo do — Xq 


2 


1-/3. 1 + 



(3) 



(4) 



Here the first terms in the square brackets stem from elec- 
tron currents to the left {x'q < xq) and the second terms 
from those to the right (a;g > Xq) of the considered elec- 
tron at Xo ■ The approximation made in the present model 
is that the velocities (3^ E^nd (3y are assumed spatially uni- 
form. 

The equations of motion then are 



dpx 
dt 

dpy 
dt 

d'y 
'dt 

dx 

It 
dy 
dt 



-Ex — PyBz, 



1 



J2 



P' 



f^x Ex Py Ey , 
Px = Px/j, 
Py = Pyh- 



From this we find 
d 



dt 

dpy 

~dt 



{l - Px) = {I - Px)Ex - Py{Ey ~ 
-{\- Px)Ey- Px{Ey^B,). 



(6) 



and the particle trajectories then follow from 

dx/dr — Px/n, (11) 
dy/dr ^ Py/n, (12) 

in parametric form with time given by t(r) = r + x{t) . 

The layer of finite thickness is divided into 100 sub- 
layers, and each sub-layer is assumed to move according 
to Eqs. ([Sl fTTjl . Finally, the density is calculated, using 
x{t, Xo) and 

n = no/ {dx/dxo). (13) 

From this we obtain spatial density distributions, as shown 
in Fig. [2 and O 



3 Comparison of model with PIC simulation 

The model is now applied to the reference case, already 
discussed in paper I [6j. It is compared with ID-PIC sim- 
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For a plane wave moving in vacuum in a;-direction, the 
laser fields satisfy the dispersion relation ujl = ck^ and 
ELyir) = Blz{t) with propagation coordinate t = t — x. 
For an electron moving along x{t), this implies dr/dt = 
1 — Px- Introducing k — ^ — px, making use of k = 7(1 — 
Px) = "f dr/dt, and recaUing Ex = Esx, Ey = ELy + Esy, 
Bz = Blz+Bsz, we derive, after some algebra, the coupled 
equations for k(t) and Pyir): 



^ - p 
dr ~ 



dPy 

dr 



Nido~xo)j^ 



9 ' 

Py 



= -E 



Nd, 



Ly 



oPy 



(7) 
(8) 
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They are solved by numerical integration. From k — "f—px, 
7^ = 1 + p'^ + Py, and the solutions Pyir) and k{t), one 
can obtain 



7(t) = 
Px{t) 



1 



2k 



pI 



2k 



K 

2' 

K 

2' 



(9) 
(10) 



Fig. 1. Density profiles of electron layer n{x) (thick lines) 
inserted into electric field ELy(x) of single-cycle laser pulse 
(thin lines) at times (a) one and (b) three laser cycles after 
first interaction with foil. Model results (solid) are compared 
with ID-PIC results (dashed). Density is normalized to critical 
density and electric field to Eq — mcLOL/e (see text). 
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ulation in Figs. [T] and [21 A planar foil with density — 
no/ricrit = 159 and eg = Nd{) = 1 is irradiated by a single- 
cycle laser pulse E^y — ao sinr (0 < t < 27r) with ap — 5. 
At this high intensity, we consider foil ionization to hap- 
pen very rapidly; actually we use the approximation of 
an initially completely ionized plasma layer. The ID-PIC 
simulations have been performed using the code LPIC++ 

The results in Fig. [1] show profiles of laser field and 
electron density at times t/ro = 1 and t/ro = 3 after 
interaction with the foil. The laser pulse has expelled all 
electrons from the foil. They are seen as a dense layer 
carried along by the first half-wave of the pulse, while foil 
ions are considered immobile and are located at x = 
(not shown). Apparently, the electron layer is transparent 
to the light, the front has penetrated the layer, but is 
depleted due to interaction with the layer. Different from 
the results in paper I, the analytical model [1 is now well 
reproducing the ID-PIC results. There are two reasons for 
this agreement: (1) correct account for the self-radiation 
of the electron layer which may be viewed as the effect of 
forward Thomson scattering reducing the wave amplitude 
in front of the layer; (2) correct account of initial electron 
acceleration inside the ion layer. 

A conspicuous deviation between model and PIC re- 
sults concerns the density spike seen in the PIC results 
on the side of the incident light pulse. It appears in sim- 
ulations with different PIC codes and may be of numeri- 
cal origin. This needs further study. Smaller deviations in 
the layer thickness develop at later times (see Fig. lb at 
t/TL = 3) and may be attributed to approximations made 
in Eqs. ^ and ([U for the self-radiation fields. 

In Figs. 2a and 2b, the temporal evolution of the elec- 
tron 7-factor is depicted. Again we find good agreement 
between model and PIC results, both for values (7) aver- 
aged over the layer (Fig. 2a) and spatial 7 gradients (Fig. 
2b). The average rises up to a maximum of (7) « 15, and 
then falls again when the layer drifts into the decelerating 
phase of the second half-wave. Notice that the maximum 
of (7) is much lower than the pure single-particle estimate 
of Jmax = 2ag — 50. This is because of the electrostatic 
and electromagnetic self- fields of the layer, well described 
by the model. 



(a) 


20 

18 
16 
14 
12 


(b) th-} 


llx=2 


t/Tj=3 , 




10 









(a) 









0.64 0.66 0.68 0.70 0.72 





0.06 


-J 


0.04 












0.02 


•5' 






O.OO 



• 


(b) 














-m- Model 




-•- PIC 





0.00 



0.01 

d A 

L 



0.02 



2 3 
t/x. 



0.68 0.70 0.72 1.56 1.60 1.64 2.46 2.52 2.58 
xlX, xJX, x/A.. 



Fig. 3. (a) Density distribution at I/tl = 1 with (solid line) 
and without (dashed line) accounting for the linear increase of 
Esx inside the ion volume (compare Eq. lU). (b) Thickness of 
electron layer d at time t/rL ~ 1 plotted versus initial thick- 
ness do; model results (solid line) compared to ID-PIC results 
(dashed line). 



4 Layer compression depending on initial 
thickness 

A particularly interesting observation is a tendency for 
layer compression during the first stage of acceleration 
stage. Compression may even overcome Coulomb expan- 
sion. It originates from the initial phase of acceleration, 
when the electrons are still inside the ion layer. There the 
Ex field, felt by an electron initially at xq and given by 
Eq. IT]), rises linearly oc (x — xq) before reaching a con- 
stant Ex value when leaving the ion volume. During this 
initial phase, electrons on the left side of the layer facing 
the incident laser pulse gain more energy than those to the 
right. This leads to layer bunching, as it is demonstrated 
for time t/T^ = 1 in Fig. [3^. Here the dashed profile corre- 
sponds to using Esx(x, Xq) = N{dQ — Xq) for all x and xq, 
while the solid profile corresponds to the correct treatment 
using full Eq. H]). 

The compression depends on the initial thickness of the 
foil. It is more pronounced when starting with thicker foils. 
In Fig. [Sb, we show the thickness d at time t/r^ = 1 as a 
function of initial thickness doj keeping the areal density 
Nd^ = 1 constant. For do/Xh = 0.02 and t/rL = 1, layer 
expansion is not only reduced relative to the reference case 
with do/XL — 0.001, but is actually shghtly compressed. 
This effect is seen in both model and PIC results. This is 
good news for experiments using laser pulses with more 
realistic shapes. In case of Gaussian shapes rather than the 
sharp-front fiat-top pulses considered here, one expects 
some extent of foil expansion before the pulse maximum 
hits. The compression then leads to high-density layers 
during the first laser cycles of layer acceleration. 



5 Expansion of a freely propagating layer 



Fig. 2. (a) Temporal evolution of the < 7 >-factor averaged 
over electron layer, (b) Spatial distribution of 7-factor over 
electron layer at different times. Model results (solid lines) are 
compared with ID-PIC results (dashed lines). 



Finally we discuss a simple solution which is important 
for the late stage of a freely propagation layer when it is 
not in contact any more with the laser field. The solution 
describes the longitudinal decay of the relativistic mirror 
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due to Coulomb explosion. Initially the layer has uniform 
density hq and plasma frequency tOp. It consists of elec- 
trons only (no ion layer) and propagates in x-direction 
with momentum 70 /3o- The equation of motion [5] then re- 
duces to 

dp^/dt = Nxo (14) 

keeping definitions and normalization the same as before. 
The coordinate xq now denotes the initial position of a 
test electron in the plane layer extending from —do to 
+do. The longitudinal momentum is 



Px{t, xo) = Polo + Nxot, 
and the corresponding energy 
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9« ^0 i f Xo t 
do To \ do To 



(15) 



(16) 



involves a characteristic time Tq = ^o/iNdo). In dimen- 
sional units, it is given by 



To 



70 



'Jpdo/c 



(17) 



Apparently, it is the relativistic plasma frequency which 
sets the longitudinal decay time. The electron trajectories 
x{t,xo) can be easily inferred from energy conservation 
7 - 70 = Ndoix -Xo), giving 



x{t, Xo) = Xo + (7 - 7o)/(iVrfo)- 



(18) 



From this the evolution of the density profile n{x,t) is 
obtained with the help of Eq. (fT3|) . Profiles of an initially 
uniform electron layer are plotted in Fig. [H In this case 
the characteristic decay time is To — 4.3 fs. 
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Fig. 4. Longitudinal Coulomb expansion of electron layer with 
initial density no — 2.7 x 10^^ cm~'^ and thickness do — 0.8 
nm, moving at 70 = 10. Density profiles are plotted for times 
t/To = (solid), 0.8 (dashed), 1.6 (dotted), and 3.2 (dash- 
dotted). The characteristic decay time is To — 4.3 fs. 
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6 Conclusion 

In conclusion, we have considered electron acceleration 
from ultra-thin foils by high-contrast ultra-short laser pulses 
in the regime in which all electrons are expelled from the 
foil. They form a dense electron layer that can be used 
as a relativistic mirror for Thomson backscattering, at 
least over a short time interval of a few laser periods after 
foil interaction. This paper supplements a companion pa- 
per discussing the backscattered spectra 6 . Here we have 
shown that the analytical model of Kulagin et al. jT] well 
reproduces corresponding particle- in-cell simulations. An- 
alytical theory describing layer dynamics is important for 
basic understanding as well as developing and analyzing 
future experiments. 

As a particular result, we have identified a regime of 
layer compression related to initial acceleration inside the 
ion volume. It is found that foils somewhat expanded ini- 
tially are superior to foils of same areal density, but thin- 
ner and denser. Also an analytical formula is given for 
longitudinal layer expansion after it disconnects from the 
driving laser field and is freely propagating. 
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